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SUMMARY 


An analytical investigation was carried out to determine the ability 
of typical solar collectors to serve as nighttime heat radiators. Such 
use might be made of the collectors during the summertime where they then 
serve as the condenser in a heat pump cycle with electrically-driven 
mechanical compressor. In winter the collectors form a solar-heated 
evaporator in the heat pump cycle. It was found that if one opens the 
upper and lower ends of a collector, large free convection currents may 
be set up between the collector surface and the cover glass(es) which can 
result in appreciable heat rejection. If the collector is so designed 
that both plate surfaces are exposed to convection currents when the upper 
and lower ends of the collector enclosure are opened, the heat rejection 
rate is 300 watts/m^ when the plate is 13°C above ambient. This is 
sufficient to permit a collector array designed to provide 100^ of the 
heating needs of a home to reject the accumulated daily air conditioning 
load during the course of a summer night. This also permits the overall 
energy requirement for cooling to be reduced by at least \ 5 % and shift 
the load on the utility entirely to the nighttime hours. 



I NTRODUCT I ON 


Using solar collectors as heat sources in heat pump cycles is a 
concept that has received favorable consideration from several designers 
of solar energy utilization systems. The corollary concept — using the 
solar collectors as nighttime radiations in heat pump cooling cycles 
is scarcely mentioned. It is evident, however, that some beneficial 
summertime use should be made of the solar collectors in such heating 
systems in order to help justify their large capital cost. Three uses 
come quickly to mind: a heat source for an absorption chiller; a source 
of domestic hot water; a heat source in a Rankine cycle power unit. 

The first use requires an extensive summer-w i nter changeover to 
utilize the collectors efficiently in the winter. If one is willing to 
sacrifice efficiency however, the collectors can be used for direct heat- 
ing (non-heat pump) in the winter with very little additional capital 
cost over that for an absorption chiller. Such systems have been built 
and tested at several locations. 

The second use Is an obvious one, but it Is readily shown that the 
collector area needed for winter heating is much greater than that needed 
for summertime domestic hot water production. This is therefore not an 
effective use of the collectors. The Rankine cycle power unit could employ 
the collectors in a system to supply part of the building’s electrical 
power needs in the summer, or it could drive a mechanical chiller instead. 

Systems have been built and operated successfully for both appli- 
cations. With present I y-ava i I ab I e technology the capital costs of a 
Rankine cycle power unit are quite high because the expander particularly 
and the condensate pump to a limited extent represent custom hardware that 
must be des i qned and fabricated for specific conditions in specific 



installations. Addition of a Rankine cycle power unit would add about 
1Q$ to the cost of the solar-assisted heat pump. 

For these reasons a less capital-intensive use for the collectors In 
the summertime was sought. It was recognized that such a use would 
probably not save as much fossil fuel energy as a more capital-intensive 
use, but it would probably be less difficult to maintain and more suited 
to individual residential applications. 

Since the solar assisted heatpump concept really implies the con- 
current use of large scale thermal storage so as to provide heat at night, 
on cloudy days, or for periods when the mechanical compressor cannot be 
operated, it was felt that the entire cycle could easily and inexpensively 
be shifted from the daytime heating mode in winter to the nighttime chill- 
ing mode in summer. Under these conditions the solar collectors become 
the freon condenser instead of the freon evaporator as when the system is 
in the heating mode. The advantages of such operation are: 

1. The system draws power from the utility at off-peak times (sunset- 
to-sunr i se) . 

2. The system COP is higher for nighttime operation than for daytime 
operation because condensing temperatures may be 6-30°C lower. 

It is conservatively estimated that because of these factors the power 
cost for rejecting a given amount of heat can be reduced by at least 15^ 
compared with a conventional air conditioning system without large-scale 
thermal storage. This estimate is based on the assumption that the solar 
collector area needed for winter heating will in fact be sufficient to 
reject all the heat necessary for summer chilling while remaining only 
slightly above ambient temperature. 

It was to investigate the validity of that assumption that the 
following analytical study was carried out. 
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COOLING MODE ANALYSIS 


An investigation into the heat transfer character i sti cs of four 
collector geometries was carried out in an effort to establish probable 
performance for solar collectors operating as nighttime radiators. The 
collector geometries considered were: 

(1) Flat plate with cover system removed (No Cover). 

(2) Flat plate with one cover (Closed Cover), 

(3) Flat plate with one cover and an opening between the plate 
and the cover at the top and bottom of the col lector 
(Open End 1 ) . 

(4) Flat plate with one cover and an opening between the plate 
and cover and between the plate and back insulation at the 
top and bottom of the collector (Open End 2). 

Selection of the collector geometries stated above stems from a 
desire to utilize a collector that will perform well both as a heat 
absorber and as a heat rejecter without requiring extensive seasonal 
modifications. Complete removal of the cover system during the summer 
is the most desirable solution from a heat transfer standpoint, but 
was deemed impractical on any but very small installations. One might 
suggest the use of a jalousie window cover system which could be opened 
during the summer and closed during the winter months. However, such an 
arrangement probably increases costs to the point of overcoming the 
collector cost advantages afforded by the heat pump concept. Additional! 
wind loads, mechanical complexity, and increased weight pose difficulties 
for a scheme of this type. 



Consideration of these problems associated with cover system re- 
moval suggests the need for a collector that can effectively dissipate 
heat with the cover system in place. Collector geometries (3) and (4) 
above are presented as a possible means of fulfilling this requirement 
without compromising heating mode performance or requiring undue mechani- 
cal complexity and additional construction costs. In both configurations 
the collector box sides at the top and bottom of the collector are hing- 
ed so that they can be opened to al low air flow across the col lector 
plate. Geometry (3) employs a conventional collector arrangement with 
the back insulation adjacent to the plate (see Figure t). In geometry 
(4) the collector plate is supported by standoffs to allow flow on both 
sides of the plate (see Figure 2). It was felt that by taking advan- 
tage of the expected free convection draft or "chimney effect" these 
types of collectors may possess adequate heat dissipation abilities. 

Cooling mode operation requires that the collector subsystem func- 
tion as the refrigerant condenser, rejecting heat to the atmosphere. 

This heat transfer is accomplished by means of convection to the sur- 
rounding air at ambient temperature and radiation from the collector to 
the night sky. The effectiveness of both of these modes of heat exchange 
is very difficult to predict in that they are Influenced by such highly 
variable parameters as wind speed, ambient temperature, atmospheric 
humidity, and cloudiness. However, It was felt that reasonable operat- 
ing limits could be obtained by assuming low wind speeds coupled with 
very clear and very cloudy sky conditions. In addition to these re- 
strictions, the following assumptions were made to simplify the heat 
transfer analysis: 
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FIGURE 1. Detail of Collector Geometry 3. 



FIGURE 2. Detail of Collector Geometry 4 
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(1) The system operation is steady state. That is, there are 
no time-dependent variables. 

(2) The absorber plate and the cover are at some average tem- 
perature. 

(3) Edge and bottom heat losses are negligible. 

(4) Thermal capacity effects of the collector plate and cover 
are negl igible. 

(5) There is negligible temperature drop across the cover. 

(6) The cover is opaque to longwave radiation (the cover 
blocks direct radiation from the plate to the sky - a good 
approximation for glass). 

(7) Back radiation from the night sky approximates that of a 
black-body at some effective sky temperature. 

(8) Radiation characteristics of the plate and cover are con- 
sidered independent of wavelength (gray-body assumption). 

The appropriate value for the effective sky temperature for various 
atmospheric conditions appears to be somewhat in doubt. For a fairly 
clear night sky Hackforth C2H suggest 273 K, while Kreith recommends a 
value of 227 K. For very cloudy sky conditions a value equal to the am- 
bient temperature at the ground would seem realistic. In an effort to 
establish suitable limits, these effective sky temperatures are consider- 
ed in combination with wind speeds of 0, 2, and 5 meters per second. 

Results of the analysis for the various geometries operating under 
the conditions stated above are presented graphically in terms of the 
heat dissipation rate as a function of the difference between the col- 
lector plate temperature and the ambient temperature (AT). Additionally, 
the effect of wind speed and the effective sky temperature at a AT of 
15°C is considered. Major elements of the heat transfer analysis for 
the collector geometries are presented in the following sections while a 
more detailed treatment for open-ended, collectors is given in the Appendix. 
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1. Flat Plate with Cover System Removed (No Cover) (Figure 3) 

The flat plate with no cover system is the most effective heat dis- 
sipator of the geometries considered. Additionally, it is the easiest 
to treat analytica i ly, the heat transfer being composed of direct radia- 
tion from the collector plate to the sky and convection from the plate 
to the surrounding air at ambient temperature. Thus, 

q , , . = q + q ( 1 ) 

^tot 1 ^r,p-s ^c,p-s 

where 


q^^^ ^ = total heat flux per unit area of collector 


q = radiative heat flux component 

^r,p-s ^ 


q _ = convective heat flux component 


If the radiation properties of the collector plate are assumed in- 
dependent of wavelength (gray), the radiative heat exchange between the 
sky (a black-body) and the plate may be found from the simple expression. 


q - & 
^r,p-s p 


a(T^ - t"^) 

P S 


(W/m^) 


( 2 ) 


where 


e = surface emittence of the collector plate 
P 

2 

a = Stefan-Bo I tzmann constant (W/m K) 

T^ = effective sky temperature ( K) 

Tp = average collector plate temperature ( K) 

The convective heat transfer from the plate to the surrounding air 
may be found from the following dimensional expression for the heat trans- 
fer coefficient recommended by Duffie and Beckman and by Meinel: 

h = (5.7 + 3.8 V) (3) 

c,p-s 
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where 


2 

h = convective heat transfer coefficient (W/m K) 

c,p-s 

V = wind speed (m/sec) 

Thus, the convective heat transfer component (q ) is given by 

^c,p-s ' 

q = (5.7 t 3.8 V) (T - T ) (W/m^) (4) 

^c,p-s p a 

where 

Tg = ambient temperature ( K) 

Thus, the total heat transferred per unit area for +he 

flat plate with no cover may be expressed in the following form: 

q+ + 1 = e a(T^ - T^) + (5.7 + 3.8 V) (T - T ) (W/m^) (5) 

^tot 1 p p s pa 

2. Flat Plate with One Cover (Closed Cover) (Figure 4) 

The heat loss components (neglecting edge and bottom loses) from a 
collector with one glass cover are shown in Figure 4. Heat exchange 
from the collector plate to the cover takes place through radiation and 
free convection heat transfer generated by the temperature gradient be- 
tween the parallel planes. Assuming no direct radiation from the col- 
lector plate to the sky, the net heat transfer from the closed cover col- 
lector ^q^Q^ 2^ heat dissipated from the cover to the night sky 

and surrounding air, the characteristics of which are as treated in the 
previous section. Thus, 

q , , „ = e q(T^ - T^) + (5.7 t 3.8 V) (T - T ) (W/m^) (6) 

TOT Zees C 0 

where 

e = surface emittance of the cover 
c 

T^ = average cover temperature ( K) 



The unknown cover temperature (T^) may be found by performing a 
heat balance on the cover as follows: 

*^c,p-c '^r,p-c ~ "^tot 2 


( 7 ) 


where 


q = convective heat transfer from the collector plate to 

^c,p-c 2 

the cover (W/m ) 


q = radiative heat transfer from the collector plate to 

^r,p-c 2 

the cover (W/m ) 


The radiation component for gray-body parallel planes is given by: 


a(T^ - t4) 

P c 


_L + _L _ 1 

£ 

c p 


The convective component is given by: 


q = h (T - T ) 

^c,p-c c,p-c p c 


( 8 ) 


(9) 


where 

2 

h = convection heat transfer coefficient (W/m K) 

c,p-c 

Determination of the convective heat transfer coefficient between 
the plate and cover has been the object of several studies, the most 
cited of which is by Tabor, presented in 1958. However, a more recent 
work by Hollands, et al. presents a correlation for calculating Nusseit 
numbers for inclined, enclosed air spaces of high aspect ratio, having 

5 

0 <_ 6 £ 60 degrees and 0 < Ra £ 10 , where 0 is the ti It angle and Ra is 


the Rayleigh number. This correlation is a closed form expression which 
has the obvious advantage of suitability for the digital calculations 
used in this analysis. The equation is as follows: 
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Nu = 1 + 1.44 


1 - 


1708 


Ra cos 0 


where 


Nu = 


h d 

P.^.P. S 


1 - 


(sin 1.8 1708 


Ra cos 
0 _< 60 degrees 


' Ra cos 0 
5830 


- 1 


(10) 


( 11 ) 


(3 ^ P 


ky J 



Ra = I ^1 1 (T^ - T ) d- 

P ^ 


d = distance between the plate and cover 

k = thermal conductivity of the air 

g = acceleration due to gravity 

p = dens i ty of air 

C = specific heat of air 
P ^ 

8 = coefficient of thermal expansion of air 


( 12 ) 


and 


y = dynamic viscosity of air 

Thus, with the use of the expressions given above, the cover tem- 
perature (T^) may be determined through an iterative process for given 
values of the collector plate temperature ^T^), the ambient temperature 
(Tg), and the effective sky temperature (T^)- The total heat dissipa- 
tion rate for the closed cover col lector may then be found from equa- 
tion (6) . 


3. Flat Plate with One Cover and an Opening Between the Plate and Cover 
at the Top and Bottom of the Collector (Open End 1 - Figure 5) 

In addition to the radiative and convective heat transfer considered 
NOTE; Starred (*) brackets go to zero when negative. 
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in the previous section, the analysis of this collector geometry is com- 
plicated by the natural convection flow through the air gap between the 
plate and cover producing a draft or "chimney effect" as shown in Figure 
5. Air at ambient temperature enters the gap, is heated by the warm 
collector plate, which produces a buoyant force, thus causing the air to 
rise until it exits at the top. 

If it is assumed that the flow is laminar over the entire length of 
the plate and that the gap is small compared to the plate length, then 
it seems reasonable to assume that the velocity and temperature profiles 
are fully developed at the exit. Additionally, if the profiles are con- 
sidered to be essentially constant the momentum and energy 
equations may be solved exactly for the velocity and temperature pro- 
files at the exit. 

With a knowledge of these relationships, a mass flow rate and a 
bulk air temperature at the exit may be calculated from which the energy 
absorbed by the air may be found. Letting "x" be the direction along 
the plate, the governing equations for two dimensional, incompressible 
flow are as fol lows: 

a) Continuity Equation 


b) X-Momentum Equation 


9u , „ 3u 


9P 


9 u 


9x - 2 

9y 


(14) 


and c) Energy Equation 


u 


IL 

9x 


+ V 


9T 

9y 



(15) 
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where 


X 

y 

u 

V 

P9s 


3x 


= direction along collector plate 
= direction across the air gap 
= velocity in the x direction 
= velocity in the y direction 
= body force in x direction 

= ai r density^ 

= pressure gradient in the x direction 
= acceleration of gravity in the x direction 
= viscosity (assumed constant) 


and 


a 


= molecular thermal diffusivlty 
Solution of these equations subject to the stated assumptions (see 
Appendix yields the following results: 

Velocity Prof i le 


_ q s I n 9 -T- ^ 

u- = ^ ^ — (T - T ) 

1 V T CO 

a 


6 6d 


+ (T 


T ) 
a 


lA _ y? 
2 2 


(16) 


where 

V = kinematic viscosity 
Temperature Profile 


Ti=^[T - T]+T 
g1 d c p-* p 


(17) 


The density (p) in the body force term is considered to vary with tem- 
perature while the density appearing elsewhere is assumed to be constant. 
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Bulk Fluid Temperature at Exit 


T 


ti^CT + T)^-T(T + T)--^T t 
15 p c a p c 15 p c 


01 


T - 2 T + T 
p a c 


(18) 


Mass Flow Through Air Gap (m^ ) 


P W g s I n 9 
S R Tqi 24 V 


(T - 2 T + T ) 
p a c 


(19) 


where 


= atmospheric pressure 

W = col lector wi dth 
R = gas constant for air 

The net convective heat transfer from the collector plate to the 

air flowing through the gap (Q ,) is given by: 

gap I 


0 - = (mJ(C )(T - T ) (W) 

gap 1 1 p o a 


(20) 


and 


, = ^ (m,)(C )(T - T ) (W/m^)(unit area)(21) 

^gap 1 A 1 p o a 


where 


A = CO I I ector area 
c 


The total heat loss 3 ) from the open ended 1 collector is 

given by: 


■^tot 3 


= q 1 + q + q 

^gap 1 ^r,c-s ^c,c-s 


(W/m^) 


(22) 


where the heat transfer components from the cover to the sky are as given 
for the closed cover geometry. Again the solution is dependent on a 
knowledge of the cover temperature (T^) which may be found from the ap- 
propriate energy balance as before. That is, 



( 23 ) 


q + q = q + q 

^r,p-c ^c,p-c ^r,c-s ^c,c-s 


where 


q = radiation from plate to cover 

r,p-c 


oCt'^ - T^) 
B 9_ 

J- + -L. 1 


(W/m^) 


(24) 


q = radiation from cover to sky 

^r,c-s ' 


4 4 

= e a(T^ - T^) 

CCS 


(W/m^) 


q = convection from cover to surrounding air 

^c,c-s ^ 


= (5.7 + 3.8 V) (T - T ) 

c a 


(W/m^) 


(25) 


(26) 


and 


q = convection from plate to cover 

^c,p-c ^ 

Determination of the convective heat transfer from the collector 
plate to the cover poses a problem due to the difficulty of obtaining a 
solution of the governing equations all along the collector length. How- 
ever the flow is assumed fully developed at the exit (v = 0) and in this 
region heat is transferred from the collector plate to the cover by pure 
conduction through the air gap. 

Since there is no knowledge of the entrance region flow character- 
istics, it is assumed that the heat convected from the plate to the cover 
is of the same order of magnitude all along the length of the collector 
as it is at the exit. This simplification yields the following result: 


q = ^ (T - T ) 

^c,p-c dp c 


(W/m^) 


(27) 


The cover temperature (T^) must again be found through an iterative 

process using the expressions given above and the total heat dissipation 
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rate for the open ended (1) collector geometry may then be found from 
equation (22). 


4. Flat Plate with One Cover and an Opening Between the Collector Plate 
and Cover and Between the Plate and Back Insulation at the Top and 
Bottom of the Col lector (Open End 2 - Figure 6) 

Analysis of the open-ended 2 collector scheme Is essentially as 

given for open end 1 collector with the addition of the free convection 

flow on the bottom side of the collector plate. That is, 

*^tot 4 '^tot 3 *^gap 2 (28) 

where 


4 = total heat transfer rate for geometry (4) 


and 


q „ = convective heat dissipated through the air gap 
below the collector plate 

Fully developed flow at the exit Is again assumed. Letting "x" be 
the direction along the plate and "y" across the air gap, the governing 
equations are identical to equations (13), (14), and (15). 

Solution to these equations subject to the appropriate boundary 
conditions yields the following results 

T 


u ^ g sin e f 1 
2 2 V T J 


(y - d, y) 


(29) 


T o = T 

g2 p 


Tno = T 

02 p 


(30) 

(31) 


P d 2 W g s I n 0 


12 R V 


T 


( 32 ) 




where subscript 2 refers to the lower air gap. Thus, 9ggp 2 


becomes 
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FIGURE 6, Heat Transfer Components in Open-End 2 Collector. 


( 33 ) 


. d^PC Wgsin9T /-T 
_ J_ 2 p ^ a P _ 1 

^gap 2 A 12 R V T T 

3 f' c p a 

and the total heat transferred for geometry (4) can be found from equa- 
tion (28). 
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RESULTS OF COOLING MODE ANALYSIS 


Results of the heat transfer analyses explained In the previous 
sections are shown in Figures 7 through 15. Property values and col- 
lector dimensions used in the analysis are shown In Table I. Note that 
the AT's are generated by holding the ambient temperature constant at 
24°C and increasing the plate temperature. Curves for the open ended 
collector geometries terminate at low AT's under certain conditions due 
to a breakdown in the free convection analysis. This occurs if the cover 
temperature drops below ambient temperature to the point that the fol- 
lowing condition is satisfied: 


T - 2T + T = 0 


p a c 

At this combination of temperatures, m^ becomes zero and the flow would 
theoretically reverse direction. As noted, this condition is only satis- 
fied at relatively low AT's and Is considered out of the range of in- 
terest from a practical heat pump operation standpoint. 

As can be seen from the graphs, the no-cover col lector geometry is 
a superior heat dissipator under nearly all the conditions shown. The 
open end 2 geometry does out-perform the no-cover collector during com- 
binations of low wind speeds, high effective sky temperatures, and fairly 
high AT's. The heat rejection ability of the open end 1 and closed 
cover collector geometry is limited except at high temperature differ- 
ences . 

The effect of wind speed and the effective sky temperature on the 
heat dissipation rate at a AT of 15°C is shown in Figures 16 and 17, re- 


spectively. Note that only the no-cover geometry is significantly sen- 
sitive to either variable. However, since the effect of wind direction 
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TABLE I. Property Values and Dimensions Used in Cooling Mode 
Analysis Calculations. 

(a) Collector Length = 2.5 meters 

(b) Collector Width (W) = 1.0 meters 

(c) Gap Distance (d) = 2.5 centimeters 

(d) Ti It Angle (6) = 45° 

5 3 

(e) Atmospheric Pressure (P^) = 1.013 x 10 Kg/m 

(f) Stefan-Bo I tzmann Constant (a) = 5.6697 x I0~^ W/m^ K 

2 

(g) Gravitational Acceleration (g) = 9.807 m/sec 

(h) Thermal Conductivity (k) = .0257 W/m K 

(i) Specific Heat (C ) = 1.005 x 10^ — - — 

P Kg K 

(i) Collector Plate Emittance (e ) = 0.98 

P 

( k) Cover Emittance (c^) = 0.94 

-5 2 

(l) Kinematic Viscosity (v) = 1.55 x 10 m /sec 

(m) Gas Constant for Air (R) = 287 — - — 

Kg K 

(n) Ambient Temperature (T^) = 24°C 



HEAT DISSIPATION RATE 






HEAT DISSIPATION RATE W/m* 


NO-COVER 


EFFECTIVE SKY TEMPERATURE = 297 K 



(COLLECTOR PLATE TEMPERATURE — AMBIENT TEMPERATURE) ®C 

FIGURE 8. Nocturnal Heat Dissipating Ability as a Function 
of the Difference Between Collector Plate Tem- 
perature and Ambient Temperature. 
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FIGURE 11. Nocturnal Heat Dissipating Ability as a Function 
of the Difference Between Collector Plate Tem- 
perature and Ambient Temperature. 
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(COLLECTOR PLATE TEMPERATURE - AMBIENT TEMPERATURE) ®C 


FIGURE 14, Nocturnal Heat Dissipating Ability as a Function 
of the Difference Between Collector Plate Tem- 
perature and Ambient Temperature. 
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HEAT DISSIPATION RATE 





HEAT DISSIPATION RATE W/m^ 


miHi iiiiM 



FIGURE 16. Nocturnal Heat Dissipating Ability as a Function of 
Effective Sky Temperature. 
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was not Included in the analysis, the results shown above should be re- 
garded with care. The free convection heat transfer component In the 
open-ended collector geometries would probably be enhanced by wind com- 
ing from most directions, although this relationship is not clearly 
estab I ished. 

The selection of a particular collector geometry for use in a 

specific installation Is Influenced by the heating as well as cooling 
mode collector performance. Of course, the desired outcome is to utilize 
the collector system that will provide a high temperature heat source in 
winter and a low temperature heat sink in summer. Consideration must 
also be given to the demand requirements on the heat pump Itself and the 
relative magnitude of loads as determined by the geographical location. 
For example, a structure in the North would have a relatively greater 
heating load and smaller cooling load than its counterpart in more 
southern regions. Additionally, heating mode collector efficiency is 
generally poorer in the North. Thus, an efficient winter absorber that 
may not necessarily be a good summer dissipator may be more cost effect- 
ive in northern locations. In the South, on the other hand, a collector 
with no cover may provide adequate winter performance, while optimizing 
summer heat dissipating ability. These general conclusions are obvious- 
ly affected by yearly heat pump operating costs versus initial collector 
cost considerations, and the choice of the most suitable option becomes 
quite involved and dependent on the characteristics of the specific ap- 
p I i cation . 

Heating mode performance for the no-cover and closed-cover collec- 
tors is shown In Figure 18 for a variety of operating conditions. The 
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HEATING MODE COLLECTOR LOSSES W/m^ 




collector heat loss Is plotted as a function of the difference between 
the collector plate temperature and the ambient temperature for various 
wind speeds. Note that for this case the AT Is generated by holding the 
collector plate temperature constant while the ambient temperature is 
decreased. A collector plate temperature of 25°C is considered in com- 
bination with wind speeds of 0, 2, 5, 10, 15, and 20 meters per second. 
Details of the heat transfer analysis for the heating mode operation are 
identical to the cooling mode analysis with the exception of the value 
of the effective sky temperature. Duffle and Beckman^ suggest the use 
of the following relation given by Swinbeck (1963) for clear days: 

T = 0.0552 T ^ 
s a 

where T and T are in degrees Kelvin. 

With the aid of the cooling and heating mode col lector performance 
curves shown above, it is felt that the selection of the collector geo- 
metry most suitable for a specific application becomes clearer. 

A brief outline of a simplified procedure for determination of the 
relative performance of the various collector options Is given below. A 
mean daily solar insolation value for January is taken from U. S. Weather 
Bureau climatological data for Greensboro, North Carolina. It is as- 
sumed that the cooling load is equal to the heating load in one case and 
then equal to one-half the heating load in the second. The general ap- 
proach taken is to size the collector area for the heating load and then 
determine the AT required to satisfy the heat dissipation rate necessary 
to meet the cooling load. Details of the analysis follow: 

^See Reference 1, pg. 76. 
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Assumptions 

2 

(a) Solar insolation on a horizontal surface = 2.32 KWH/m -day 

(b) Solar insolation on a surface tilted 45° = 3.55 KWH/m^-day 

(c) Heating load = 322 KWH/day 

(d) Can col lect heat over seven hour period 

(e) Wind speed for heating and cooling modes = 2 m/sec 

(f) Effective sky temperature for the cooling mode = 273 K 

(g) The solar collector is operated in the heating mode at a 
AT of 15°C. 

(h) Can dissipate heat over eight hour period 

Ca I cu I at ions 

(a) Heating load 

First determine the average solar insolation for the seven 
hour heat collecting period. Thus, 

2 

Average solar heat available = 3.55/7 = 507 W/m 

The heat loss under the conditions stated above can be found 

from Figure 41 for the no-cover and closed-cover geometries. 

2 

No cover heat loss = 335 W/m 

2 

Closed cover heat loss = 110 W/m 

The heat collection ability of the two geometries is then found as 
fol lows: 

No cover = 507 - 355 = 172 W/m^ 

Closed cover = 507 - 110 = 397 W/m^ 

Thus, the total collector area required is 

No cover area = 32,200/( 172) (7) = 267 m^ 

Closed cover area = 32,200/(397) (7) = 116 m^ 
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(b) Cool ing Load 

(i) Case I - Cooling load = heating load = 322 KWH/day 
For eight hour night must dissipate 332/8 = 40.2 KWH/hr 

No-cover average heat dissipation rate = 

40,200/267 = 151 W/m^ 

Closed-cover average heat dissipation rate = 

40,200/116 = 347 W/m^ 

From Figure 34 the AT's required for the heat dissipation rates 
above are as follows: 

No- cover - 1°C AT 
Closed-cover - off the graph 
Closed-cover to open-end 1 = 25°C AT 
Closed-cover to open-end 2 - 14.5°C AT 
Note that is is assumed the no-cover collector remains as is, while 
the closed cover collector may be changed to either of the other 
geometries for cooling mode oepration. 

(ii) Case II - Cooling load = one-half heating load = 

161 KWH/day 

The average heat dissipation rates are as follows: 

2 

No-cover col lector = 75 W/m 

2 

Closed-cover collector = 173 W/m 
Again, from Figure 34 the AT*s are as follows: 

No-cover - radiation alone is sufficient at 0°C AT 
Closed-cover - 23°C AT 
Closed-cover to open-end 1 - 16°C AT 
Closed-cover to open-end 2 - 10°C AT 
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For optimum performance of the heat pump a small AT for both the 
heating and cooling modes Is required. Both these objectives may be 
realized by increasing the collector area. However, this design option 
obviously leads to higher initial costs and the trade-off becomes one of 
increased initial cost versus higher operating costs. The matter is 
further complicated by the various collector alternatives, as can be 
seen even from the simplified analysis shown above. 

While the selection of the most cost - effective collector alter- 
native is a complex matter, the results of the preceding analysis do sug- 
gest the following general conclusions: 

(1) It seems probable that for many applications the open-end 
2 geometry presents a feasible means of utilizing the 
collector system as a heat sink without sacrificing heat- 
ing mode performance or substantially adding to Initial 
costs . 

(2) Due to the severe heating mode losses suffered at even 
moderate wind speeds. It is felt that the no-cover col- 
lector would be practical only in applications where the 
cooling load exceeds the heating load and where high 
levels of solar insolation are available. This alterna- 
tive would also seem the better choice if a lower heat- 
ing mode coefficient of performance is acceptable. 

(3) The open-end 1 and closed-cover collector geometries are 
unsuitable for use as heat dissipators. 

In closing, it should be recognized that the results derived above 
are approximate and their principal importance is in presenting a per- 
formance comparison of the collector geometries considered. Absolute 
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values for the heat dissipation rates of the collectors under actual 
operating conditions must obviously be derived by experimental means. 
However, the conclusions drawn above do indicate that utilization of the 
collector system on a year-round basis is a practical means of justify- 
ing the high initial costs of a solar energy installation, particularly 
if the heat pump concept is employed. 
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CONCLUSION 

The results of the present cooling mode analysis coupled with 
previously reported heat mode performance suggest that the solar assisted 
heat pump system is a viable alternative to the more conventional means 
of solar energy utilization. Thus, even though economic implications 
have not been fully explored, it is felt that construction of a larger 
scale prototype is certainly justified. The concluding remarks given 
here will therefore be limited to suggestions with regard to design con- 
siderations of such a system. 

Since the response of the refrigerant flow control system is con- 
sidered the most critical element in optimizing the performance of a 
solar assisted heat pump, it is recommended that externally driven 
compressors be specified for installation in a multiple compressor 
arrangement. This would allow maximum flexibility of control of the 
compressor displacement without requiring extensive modifications. 
Furthermore, the fine adjustment capability thus achieved would permit 
the establishment of optimum flow rates for various operating conditions 
and even allow for refrigerant substitutions. Additionally, most 
externally driven compressors are provided with a sight glass in the 
pump for monitoring the oil level, an important consi deration in view 
of the difficulties encountered in the NCSU system. The use of com- 
pressors of this type would also have the benefit of the relative ease 
with which periodic inspections of the valves and piston rings could 
be made. 
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In conclusion, it is hoped that the material presented in this 
report wi I i stimuiate sufficient interest to encourage further research 
and development of solar assisted heat pumps. As noted earlier, the 
NCSU system was built on a limited budget and was undertaken primarily 
as a design exercise for undergraduate engineering students. Therefore, 
significant performance improvement is considered possible. 
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APPENDIX 


OPEN-ENDED COLLECTORS HEAT TRANSFER ANALYSIS 


A. Open-End 1 Col lector 

For sake of clarity, the heat transfer components involved in the 

open-end 1 collector geometry are shown in Figure 19. Treatment of the 

radiation components, q and q , and the convection components, 

r ^ p "" c r* ^ c* s 

q and q _ > Is given in the cooling analysis. Therefore, con- 

C^p^C C^C"S 

sideration here is limited to the development of the natural convection 
draft problem. 

As mentioned in the cooling mode analysis, air at ambient tempera- 
ture (T ) enters the opening at the bottom of the collector, is heated 
by the warm plate and exits through the top opening. The net energy 
transferred from the plate to the air (Q .) is a function of the air 
mass flow rate and the bulk fluid temperature^ exit. That 

is. 


Q 


gap 




It is assumed that the flow is essentially 
developed laminar flow exists at the exit. 

u = u(y) only 


- T ) (A.l) 

a 

incompressible and that fully 
Thus, 


and 


T = T(y ) only 

The bulk fluid temperature is defined as: 


Also known as mixed mean fluid temperature. 
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(A. 2) 


T 


01 


u(y) T(y) dy 


u(y) dy 


and the mass flow is defined in terms of the mean fluid velocity u as 


m^ = p A u 
1 gap 


(A. 3) 


where 


1 

^ = d 


u(y) dy 


(A. 4) 


and 


A = dW 

gap 


(A. 5) 

As can be seen from the relations shown above, determination of 


0 1 requires the solution of the governing flow equations for the tern- 

gap 1 

perature and velocity profiles. The governing equations are as follow: 
(a) Continuity Equation 


3x 3y ^ 


(A. 6 ) 


(b) X-Momentum Equation 


3u , ,, 9u 


(c) Energy Equation 


- ^ n 4. ,, 


(A. 7) 


3T ^ 3T 3^ 

u v" + V - 5 — = a — ~ 

3x 3y 9 y 2 


(A. 8 ) 


Since the flow is generated by thermal effects, the equations are 
coupled and the solution for the velocity profile requires knowledge of 
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the temperature distribution. Consider the energy equation. If it as- 
sumed that the temperature profile is essentially constant In the region 
2 

near the exit the analysis is greatly simplified and the energy equa- 
tion reduces to an ordinary differential equation. The boundary condi- 
tions are constant wall temperatures, Tp at the plate and at the cov- 
er. Thus, 3T/9x “0, v = 0 for fully developed flow and constant pro- 
files and 

T = T @ y = 0 
P 

T = T @ y = d 
c ^ 


Therefore, 


i ntegrat i ng 


and 


but 


and 


Thus, 



9y 


T 



+ C 


2 


T = T 


P 


y = 0 


T = T @ y = d 
c ' 




- T ) + T 
P 


P 


(A. 9) 


(A. 10) 


(A.11) 


(A. 12) 


2 


See Reference 


3, pp. 114-117. 
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Note the resulting temperature profile is linear in y. Next con- 
sider the continuity and x-momentum equations. Since the flow Is fully 
developed (v = 0), the continuity equation yields the following result: 
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3x 


= 0 


(A 


Thus, the x-momentum equation becomes 


. 3P ^ 3^u 


(A 


where - p g = -the body force in the x-di recti on. 

X 

The pressure force term under these conditions is the external ' 
drostatic" pressure gradient due to changes in elevation. Thus, 


3P _ 

3x *^a ^x 


.A 


Substitution into equation (2.12) yields 


0 = g^Cpg - p) + y 


3^u 

3y^ 


(A, 


The density difference (p - p) may be approximated in terms of 

9 

volume coefficient of expansion, as follows: 




'3V' 


.9L 


. V - V 
J a 

V T - T 
a a 


(A 


or 


B = 


P - P 
p(T - T ) 

a 


(A. 


Thus, 


(Pg - p) = p(T - Tg)3 


(A 


Substitution of equation (A. 17) into equation (A. 14) and dividing by 
yields. 


.13) 

.14) 

”hy- 

.15) 

,16) 

the 

17) 

18) 

19) 

P 



(A. 20) 


.2 

0 = S(T - T ) + v-^ 

^ ^ 3y^ 

But, g^ = g sin 0 for the geometry given and T = ■^ (T^ - T^) from the 
equation (A. 10). Thus, after rearranging, the x-momentum equation be- 
comes 


3 u _ _ B q si n 6 [y_ 
.2 V 1 6 

3y ‘- 


(T - T ) + T - T 


(A. 21) 


The no slip boundary conditions are 

u = 0 @ y = 0 
u = 0 @ y = d 

Integrating equation (2.19) twice yields 
B q s i n 9 


(A. 22) 


ui = - 


V 


3 2 

^(T - T)+^(T - T) 

6d c p 2 p a 


+ (A. 23) 


Application of the boundary conditions results in the following 

(A. 24) 


r = B q sin 0 
^1 V 


-j(T - T)+4(T - T) 

,6 c p 2 p a 


and 


C 2 = 0 


So the velocity profile becomes 

= B g sin, e ^ ( _ T ) fid _ yI' 

^1 V c p [6 6dJ 


+ (T„ - T ) 
P a 


'll _ if ] } 
[2 2}^ 


(A. 25) 


If the air is considered an ideal gas the volume coefficient of expansion 
may be expressed as fol lows 

P. 




1 


(T - T ) 
a 


from equation (A. 16) and 
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p 

p = ^ 


from the ideal gas law. So 


or 






T_ 


1 ) 


a 


(A. 26) 


With a knowledge of the approximate velocity and temperature pro- 
files derived above the solution for the bulk fluid temperature and the 
mass flow rate may be carried out by substituting the appropriate func- 
tions in equations (A. 2), (A. 3), and (A. 4). Performing the necessary 
integrations results in the following expressions 


T 


01 


q^(T +T)^-T(T + T)-^T T 

15 p c a p c 15 p c 

T - 2 T + T 
p a c 


(A. 27) 


- , d Bg sin e „ . 

1 24 V P a c 


(A. 28) 


and 


p Wd^BgsinS 

m. = -2 (T - 2 T + T ) (A. 29) 

1 24 V Tg p a c 

where is the density evaluated at the bulk fluid temperature which 

can be expressed with the use of the perfect gas law as follows 


P 


o 



(A. 30) 


where P is the atmospheric pressure, assumed constant. Thus, after sub- 
stituting for 3 from equation (A. 26), the flow rate, becomes 
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(A. 31) 


P W d g s i n 0 
• _a 

"’l R 24 V T. 


(T - 2 T + T ) 
p a c 


01 “ ■ ' a 

and the net energy absorbed by the air in flowing through the collector 
is given by equation (A.1). For a collector of characteristic dimen- 
sions, the average heat transferred per unit area may be found by divid- 
ing by the collector area. Thus 

1 


. (m. )(C )(T.. - T ) 

A^ 1 p 01 a 


(A. 32) 


*^gap 1 

Note that the subscript 1 refers to the gap above the collector plate. 


B. Open-End 2 Collector 

As in the previous section, the discussion here is limited to the 
convection problem through the air gap. The principal heat transfer com- 
ponents are shown in Figure 20. Since It is assumed that the collector 
plate temperature is held constant, the heat transfer components above 
the plate are unaffected by the convection through the lower air gap. 
Therefore, consideration here is given only to the development of the 
free convection below the collector plate. Again it is assumed the flow 
is fully developed, laminar and essentially incompressible in the region 
near the exit. The constant profile assumption is also required for 
simplification. Thus, the governing equations reduce to the following 

(a) Continuity 

|^=.0 (A.33) 

(b) X-Momentum 

3y^ ® 
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FIGURE 20. Heat Transfer Components in Open-End 2 Collector. 



(c) Energy Equation 


(A. 35) 


Note the equations are coupled as before, and the solution first requires 
determination of the temperature profiie. The appropriate boundary con- 


ditions are as foi lows 


T = T @ y = 0 
P 


0 @ y = 


^1y 


C, = 0 , = Tp 


T o = T 

g2 p 


(A. 36) 


(A. 37) 


Note the temperature is constant all across the air gap and the 


bulk fluid temperature 1^2 equal to the plate temperature of 

^02 = Tp 


(A. 39) 


Substitution of equation (A. 37) into the x-momentum equation yields the 
fol lowing 


3 u _ g g sin 0 
.2 V 

3y 


[T - T ] 
•■a p-' 


(A. 40) 
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(A. 41) 


9 u _ q sin 9 

a 2 V 

9y 




T n 
-E. 

■'a. 


and after integrating 


u = 


q s i n 9 


2 V 


T 1 
a-' 


y + y + C, 


The no slip boundary conditions are 

u = 0 @ y = 0 

u = 0 @ y = d 2 

Thus 

T 


q s I n 9 

u = 

2 2 V 


’-t" 

a- 


(y - d 2 y) 


Substitution of these profiles into the expressions for u and 


_ d g sin 9 (-T 

u — — ^ — 1 

^2 12 V T ' 


^o R T 


m2 = 


P d:; W g sin 9 r 
a 2 ^ 


12 R V 


11 


• a p-- 


Thus, the net energy absorbed by the air flowing through 
gap, on a unit collector area basis, is given by 


■^gap 2 


d:; P C W g sin 9 T 
_ 2 a p ^ a 


A 12 R V 
c 


T 


p L a 


T 


(A. 42) 


(A. 43) 
m yields 

(A. 44) 
(A. 45) 

(A. 46) 
the lower 

(A. 47) 
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